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Dark matter annihilations taking place in nearby subhalos could appear as gamma-ray sources 
without detectable counterparts at other wavelengths. In this study, we consider the collection of 
unassociated gamma-ray sources reported by the Fermi Collaboration in an effort to identify the most 
promising dark matter subhalo candidates. While we identify 24 bright, high-latitude, non-variable 
sources with spectra that are consistent with being generated by the annihilations of ~ 20-70 GeV 
dark matter particles (assuming annihilations to bb) , it is not possible at this time to distinguish these 
sources from radio-faint gamma-ray pulsars. Deeper multi-wavelength observations will be essential 
to clarify the nature of these sources. It is notable that we do not find any such sources that are 
well fit by dark matter particles heavier than ~100 GeV. We also study the angular distribution of 
the gamma-rays from this set of subhalo candidates, and find that the source 3FGL J2212.5+0703 
prefers a spatially extended profile (of width ~ 0.15°) over that of a point source, with a significance 
of 4.2cr (3.6cr after trials factor). Although not yet definitive, this bright and high-latitude gamma- 
ray source is well fit as a nearby subhalo of m x ~ 20-50 GeV dark matter particles (annihilating to 
bb) and merits further multi-wavelength investigation. Based on the subhalo distribution predicted 
by numerical simulations, we derive constraints on the dark matter annihilation cross section that 
are competitive to those resulting from gamma-ray observations of dwarf spheroidal galaxies, the 
Galactic Center, and the extragalactic gamma-ray background. 

PACS numbers: 95.35.+d, 95.95.Pw, 07.85.-m, FERMILAB-PUB-15-124-A 


I. INTRODUCTION 


Numerical simulations reveal that dark matter struc¬ 
tures form hierarchically, beginning with the smallest ha¬ 
los, and gradually merging to create ever larger systems, 
including the halos that host galaxies and galaxy clus¬ 
ters p. As a consequence of this process, dark mat¬ 
ter halos are predicted to contain very large numbers 
of smaller subhalos. In the case of the Milky Way, the 
largest members of this subhalo population include the 
few dozen known dwarf galaxies and the Large and Small 
Magellanic Clouds. These systems are exceptional, how¬ 
ever, and reflect only the very small fraction of subhalos 
that were large enough to capture significant quantities 
of gas and form stars. The vast majority of the Milky 
Way’s subhalos harbor no significant quantities of bary- 
onic matter and cannot be detected by optical surveys. If 
dark matter particles annihilate at a rate similar to that 
expected of a simple thermal relic, however, the nearest 
and most massive subhalos could generate an observable 
flux of gamma-rays, appearing as a population of uniden¬ 
tified gamma-ray sources. 

The Fermi Collaboration has recently released a new 
catalog of gamma-ray sources, known as the Fermi- 
LAT Third Source Catalog, or 3FGL [2]. Along with 
many identified objects, this catalog contains 992 sources 
that have not been associated with emission observed at 
other wavelengths. These sources almost certainly in¬ 
clude many presently unidentified blazers, pulsars, and 


other astrophysical objects^ If the dark matter consists 
of annihilating particles with weak-scale masses, how¬ 
ever, we should also expect a relatively small number 
of these unassociated sources to be dark matter subha¬ 
los mm- For example, we estimate that for a 100 GeV 
dark matter particle with an annihilation cross section 
of crv ^ 2 x 10 -26 cm 3 /s, the 3FGL should contain on 
the order of ~10 sources which are, in fact, dark matter 
subhalos. 

The challenge, of course, lies not in merely detecting 
such subhalos, but in differentiating them from the much 
more numerous conventional unidentified sources. With 
this goal in mind, we expand upon previous work [5HEE] in 
a number of ways. First of all, we employ Fermi’s recently 
recently catalog of gamma-ray sources, the 3FGL [2 . We 
also make use of the current Fermi dataset to extract a 
finely binned spectrum from each unassociated source, 
allowing us to make detailed comparisons with the pre¬ 
dictions of various dark matter models. From this infor¬ 
mation, we place constraints on the dark matter annihila¬ 
tion cross section that are competitive with those derived 
from dwarf galaxies nans], the Galactic Center m, and 
the extragalactic gamma-ray background Hang, we 
also identify a collection of 24 bright and high-latitude 
gamma-ray sources with dark matter-like spectra. We 
further investigate this list of prospective dark matter 


1 In the Galactic Plane, some 3FGL sources also remain unas¬ 

sociated due to challenges in discriminating between multiple 

astrophysical objects in the same region of the sky. 
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subhalo candidates by testing their angular distribution, 
and find evidence of spatial extension from the subhalo 
candidate source 3FGL J2212.5+0703. 

The remainder of this paper is structured as follows. 
In Sec. [TTJ we describe our calculation of the distribution 
of dark matter subhalos and their predicted gamma-ray 
fluxes. In Sec. m we discuss the characteristics of the 
unassociated source population presented in the 3FGL 
catalog, and place cuts on variability and galactic lat¬ 
itude in an effort to separate prospective subhalo can¬ 
didates from astrophysical sources. In Sec. m we de¬ 
scribe our analysis of the Fermi data and the determina¬ 
tion of the spectra from the 3FGL’s unassociated sources. 
Focusing on the most promising subhalo candidates, we 
test these sources for indications of spatial extension in 
Sec.[V] In Sec. VI, we use the number of observed subhalo 
candidates to derive constraints on the dark matter anni¬ 
hilation cross section, finding limits that are competitive 
with the strongest constraints from other gamma-ray ob¬ 
servations. Finally, in Sec VII, we summarize our results 
and comment on the prospects for future study. 



FIG. 1. The number of high latitude (|6| > 20°) dark matter 
subhalos predicted by our calculations to be bright gamma- 
ray sources (F 1 > 7 x 10“ 10 cm -2 s -1 , E 1 > 1 GeV), as a 
function of the annihilation cross section, for three values of 
the mass. Here, we have assumed annihilations that proceed 
to bb. 


II. GAMMA-RAYS FROM NEARBY DARK 
MATTER SUBHALOS 

The prospects for observing gamma-rays from nearby 
dark matter subhalos depend not only on the character¬ 
istics of the dark matter particle itself, but on the local 
abundance and density profiles of the subhalos. To esti¬ 
mate these quantities, we make use of the results of the 
Aquarius Project, which has provided the highest resolu¬ 
tion simulations to date of the dark matter subhalo pop¬ 
ulations found within the halos of Milky Way-like galax¬ 
ies, identifying hundreds of thousands of subhalos, with 
masses as low as 3.24 x 1O 4 M 0 El- The mass distribu¬ 
tion of these subhalos takes the form of dN/dM oc M -1-9 , 
with an overall normalization that corresponds to 13.2% 
of the Milky Way’s total mass in dark matter. 

In our calculations, we follow the approach of Ref. 0, 
to which we direct the reader for further details. We 
will, however, repeat the main elements of our calcula¬ 
tion here. Firstly, based on the results of the Aquarius 
simulation, we adopt the following distribution for sub¬ 
halos in the local volume of the Milky Way: 


each subhalo (prior to tidal effects), defined as the ra¬ 
tio of the virial and scale radii, c = r v i r /r_ 2 , we adopt 
values as presented in Ref. [18], with subhalo-to-subhalo 
variations modeled by a log-normal distribution with a 
dispersion of a c = 0.24 m- In this respect, we depart 
from the mass-concentration relationship used in Ref. [5], 
as based on the results of Ref. [20]. This update of the 
mass-concentration relationship reduces the number of 
detectable subhalos (in gamma-rays) by a factor of ap¬ 
proximately 4.6 relative to that presented in Ref. [5]. 

The population of subhalos detectable in gamma-rays 
is dominated by the most massive and nearby members 
of this population. In our calculations, we include sub¬ 
halos with masses up to 10 7 M o . We have chosen to limit 
our calculation to subhalos below this mass because we 
expect many of the more massive subhalos to contain sig¬ 
nificant quantities of baryons (stars and/or gas) and thus 
would be identified with dwarf spheroidal galaxies. 

The differential gamma-ray spectrum per solid angle 
from dark matter annihilations within an individual sub¬ 
halo is given by: 
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By integrating this expression between M = 3.24 x 
10 4 M 0 and 1O 7 M 0 , we find a local mass density in 
subhalos of 5700 M 0 /kpc 3 (0.00022 GeV/cm 3 ), corre¬ 
sponding to approximately 0.054% of the overall local 
dark matter density (in good agreement with Fig. 12 of 
Ref. [12]). 

For each individual subhalo, we adopt an Einasto pro¬ 
file with a = 0.16, truncated by the effects of tidal strip¬ 
ping, such that only the innermost 0.5% of the mass 
remains intact @E|. For the initial concentration of 


where m x is the mass of the dark matter particle, ( crv ) 
is the annihilation cross section, and dN 1 /dE 1 is the 
gamma-ray spectrum produced per annihilation, which 
we calculate using PYTHIA 8 [23. The integral of the 
density squared is performed over the line-of-sight, D 
is the distance to the center of the subhalo, 0 is the 
angle to the center of the subhalo, and r(#,D,/) = 
V-D 2 + « 2 -2Z)Zcos6». 

In Fig. [T| we plot the number of high latitude (\b\ > 
20°) dark matter subhalos that are predicted to be bright 
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Variability Index 



FIG. 2. Left: The distribution of the variability index observed among the 3FGL’s unassociated source population. The 
dashed curve represents the statistical distribution predicted for a population of non-variable sources (a y 2 distribution with 
47 degrees-of-freedom). Right: The latitude distribution of the 3FGL’s unassociated sources. Based on a simple disk+isotropic 
population model (dashed curve), we estimate that approximately 20% of these sources are part of an isotropic population. In 
our analysis, we will focus on the unassociated sources with a variability index less than 80 and that are located at \b\ > 20° 
(| sin b\ > 0.342). 


gamma-ray sources (F 7 > 7 x 10“ 10 cm -2 s _1 , E 7 > 1 
GeV, which is well above the threshold for inclusion in 
the 3FGL catalog), as a function of the annihilation cross 
section, for three choices of the mass (for the case of an¬ 
nihilations proceeding to 66). For an annihilation cross 
section near the value predicted for a simple thermal relic 
(<jv rsj 2 x 10 -26 cm 3 /s), one expects such bright subhalos 
to exist only for relatively low values of the dark matter 
mass, m x <100 GeV. For lighter masses, however, one 
predicts that Fermi should have already detected gamma- 
ray emission from several such subhalos, constituting a 
population of presently unidentified sources without as¬ 
sociated emission at other wavelengths. 

We emphasize that although non-negligible uncertain¬ 
ties do exist in our calculation of the number of dark 
matter subhalos above a given gamma-ray flux, the over¬ 
all conclusions of these calculations are generally robust. 
Reasonable changes in our assumptions regarding the lo¬ 
cal number density of subhalos and dark matter distri¬ 
bution within subhalos could plausibly increase or de¬ 
crease the predicted number of gamma-ray sources by up 
to a factor of a few. That being said, we have made 
a number of choices which we consider to be conserva¬ 
tive, and thus expect our calculations to reflect a fairly 
realistic, albeit somewhat low, estimate for the number 
of such sources that might be observed. For example, 
steeper density profiles (as suggested by the Via Lactea 
II simulation [22] ) could increase the predicted number 
of observable gamma-ray subhalos. Furthermore, if we 
had included subhalos larger than our cutoff of 10 7 M©, 
greater numbers of such sources would be predicted. We 
have also neglected any annihilation “boost factors” that 
might result from substructure within individual dark 
matter subhalos. 


III. UNIDENTIFIED SOURCES IN THE 3FGL 
CATALOG 

The Fermi Collaboration has recently released their 
third source catalog (3FGL), based on their first four 
years of data [2]- This catalog contains 3033 sources, 
each detected with greater than approximately 4cr signif¬ 
icance. About half of these sources have been identified 
as, or associated with, known active galaxies (including 
BL Lacs, flat spectrum radio quasars, and other classes 
of active galaxies). A smaller, but not insignificant num¬ 
ber, of these sources have also been associated with galac¬ 
tic objects, including pulsars, supernova remnants, and 
globular clusters. Of the 3033 sources contained in the 
3FGL, 992 have not yet been identified or associated with 
emission observed at other wavelengths. It is among this 
subset of 3FGL sources that we could potentially find a 
population of dark matter subhalos. 

The 3FGL provides information about each source that 
we will use to refine our search for dark matter subha¬ 
los. First, the Fermi Collaboration has tested each source 
for variability, as can be exhibited by some classes of 
astrophysical sources, but not by dark matter subhalos. 
Note that this test is performed by dividing the data into 
month-long temporal bins, and thus is not sensitive to 
variations taking place over shorter timescales, such as is 
observed from pulsars, for example. For each source, the 
Fermi Collaboration reports a value for its “variability 
index”. In the left frame of Fig. [2] we plot the distribu¬ 
tion of this quantity observed among the 3FGL’s unas¬ 
sociated source population. The dashed curve represents 
the statistical distribution predicted for a population of 
non-variable sources (a x 2 distribution with 47 degrees- 
of-freedom) [2]. For variability indices lower than ap- 
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proximately 80, these results are in good agreement, and 
thus provide no evidence of a variable population. In ad¬ 
dition, we identify 18 unassociated 3FGL sources with a 
variability index greater than 80. In our analysis, we re¬ 
move these 18 sources from our list of potential subhalo 
candidates. 

In contrast to other galactic gamma-ray sources, which 
tend to be concentrated near the Galactic Plane, bright 
dark matter subhalos are predicted to be approximately 
isotropically distributed on the sky. We can therefore use 
galactic latitude as an indicator of the likelihood that a 
given gamma-ray source is a dark matter subhalo. In the 
right frame of Fig. [2| we plot the latitude distribution of 
the 3FGL’s unassociated sources. As a histogram binned 
in sin(b), an isotropic distribution would be flat, whereas 
the observed distribution clearly includes a population of 
sources that is concentrated around the Galactic Plane. 

To estimate how this distribution breaks down into 
different spatial populations, we model the total distri¬ 
bution of sources as the sum of an isotropic component 
and a component with a thick disk-like distribution. For 
the disk-like distribution, we adopt n oc exp(— z/zq), and 
assume simply that individual sources are detectable out 
to a common distance of d max . Allowing the normaliza¬ 
tions of the isotropic and disk-like components to vary, 
along with the parameters zq and d max , we find that 
the observed distribution is best-fit by zo/^max — 0.1 
and ni S otropic/^totai — 0.45 at b = 90° (shown as a 
dashed curve in the right frame of Fig. [2|. In this fit, 
we have neglected the innermost two bins of the distri¬ 
bution, which are highly biased by the incompleteness 
of multi-wavelength AGN catalogs in this region of the 
sky El- This simple model suggests that roughly ^20% or 
^200 of Fermi’s 992 unassociated sources are part of an 
isotropic population, consisting of extragalactic sources, 
and perhaps a small number of nearby dark matter sub¬ 
halos. In our analysis, we will focus on Fermi’s unasso¬ 
ciated sources located at \b\ > 20°, allowing us to limit 
contamination from galactic astrophysical sources. 


IV. FERMI DATA ANALYSIS 

To calculate the gamma-ray spectra from Fermi’s unas¬ 
sociated sources, we utilize data taken over approxi¬ 
mately 6.4 years of Fermi-LAT observations^] utilizing 
the Pass 7 Reprocessed photons in the energy range of 
100 MeV to 100 GeV. We exclude events arriving at a 
zenith angle greater than 100°, as well as those which do 
not pass the “Source” photon data selection. We exclude 
events recorded while the instrument was not in science 
survey mode, when the instrumental rocking angle was 
>52°, or when the instrument was passing through the 
South Atlantic Anomaly. We examine the photons ob¬ 


2 MET range 239557417 - 442851954 


served within a 14° x 14° box centered around the loca¬ 
tion of each source, and divide the photons into 140x140 
angular bins and 15 evenly spaced logarithmic energy 
bins. We analyze the instrumental exposure throughout 
the region-of-interest using the P7REP_S0URCE_V15 instru¬ 
mental response functions. 

In our analysis, we employ the latest model for diffuse 
galactic gamma-ray emission, gll_iem_v05_revl. f it, 
the latest isotropic emission template for the Source pho¬ 
ton data selection iso_source_v05.txt, and include all 
3FGL sources which lie within our region-of-interest. We 
allow the normalization of each source to float indepen¬ 
dently in each energy bin, and do not impose any pa¬ 
rameterization on their spectral shape. Our model also 
includes 3FGL sources which lie nearby (but outside) of 
the region-of-interest. To calculate the best fit flux from 
each source (in each energy bin), we use the Fermi-LAT 
gtlike code, utilizing the MINUIT algorithm. We also 
calculate the upper limit of the source flux using the 
pyLikelihood UpperLimits tool in each energy bin, and 
present only the (2 a) upper limit whenever the flux is 
smaller than twice the calculated flux error. 

In Figs. [3] and [4] we plot the gamma-ray spectra from 
the 14 brightest (F 7 > 2 x 10 -9 cm -2 s -1 , E 7 > 1 GeV) 
non-variable (variability index < 80), unassociated 3FGL 
sources located outside of the Galactic Plane (\b\ > 5°). 
The 9 sources shown in Fig. [3] are reasonably well fit 
(x 2 < 25) by dark matter annihilating to 65, and the 
dashed curves represent this prediction for the best-fit 
value of the dark matter mass. The 5 sources shown in 
f^. a in contrast, are not well-fit for any choice of the 
dark matter mass. 

Of the well-fit sources shown in Fig. [3] we note that all 
9 prefer dark matter masses in the range of roughly 20 to 
70 GeV. This is further explored in Fig. [5] where we plot 
the galactic latitude and best-fit dark matter mass (for 
annihilations that proceed to bb) for each of the unas¬ 
sociated 3FGL sources with F 7 > 1.0 x 10 -9 cm -2 s -1 
(E 7 > 1 GeV), variability index < 80, and |6| > 5°. 
The error bars around each point represent the range of 
dark matter masses for which the fit to the spectrum 
yields x 2 < 25. We note two things about this plot. 
First, the distribution of these sources is not isotropic, 
and a component concentrated around the disk is clearly 
evident. Among the 12 of these sources with |6| > 20°, 
however, the distribution is consistent with isotropy. Sec¬ 
ond, nearly all of these sources favor dark matter masses 
in the range of approximately 20 to 70 GeV. While this 
could represent an indication of dark matter subhalos 
with m x ^ 30 — 50 GeV, the spectral shape in question 
is similar to that observed from many gamma-ray pul¬ 
sars. 

In Table [T] we list what we consider to be Fermi’s most 
promising dark matter subhalo candidates. This includes 
all unassociated and non-variable sources with |6| > 20°, 
<F 7 > 7 x 10 -1 ° cm -2 s -1 (> 1 GeV), and that are well- 
fit (x 2 < 25) by annihilating dark matter for at least 
one value of the dark matter’s mass (assuming annihi- 
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FIG. 3. The gamma-ray spectra from 9 of the brightest (F 1 > 2 x 10 9 cm 2 s E 1 > 1 GeV) non-variable (variability index 
< 80), unassociated 3FGL sources located outside of the Galactic Plane ( \b\ > 5°). Each of these 9 sources are reasonably well 
fit (x 2 < 25) by dark matter annihilating to bb : and the dashed curves represent this prediction for the best-fit value of the 
dark matter mass. 
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FIG. 3 (Cont.). 
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FIG. 4. As in Fig. [ 3 ] but for the 5 bright (F 1 > 2 x 10 -9 cm -2 s _1 , E 1 > 1 GeV) non-variable (variability index < 80), 
unassociated 3FGL sources located outside of the Galactic Plane (\b\ > 5°) that are not well-fit by any dark matter particle 
annihilating to bb. 

















































Source Name 

<f> 7 (cm 2 s 1 ) 

b (deg) 

m x (GeV), best-fit 

m x (GeV), Ax 2 = 4 

3FGL J0312.1-0921 

9.49xlO“ 10 

-52.36 

31.1 

12.0 - 72.3 

3FGL J0318.1+0252 

1.23xl0 -9 

-43.64 

25.5 

19.5 - 32.2 

3FGL J0456.2-6924 

7.62xl0 _1 ° 

-35.28 

10.9 

< 28.9 

3FGL J0953.7-1510 

1.25xl0“ 9 

29.61 

35.2 

26.9 - 47.8 

3FGL J1119.9-2204 

2.70xl0 -9 

36.06 

20.2 

17.2 - 23.7 

3FGL J1120.6+0713 

l.lOxlO -9 

60.69 

32.8 

23.3 - 50.5 

3FGL J1221.5-0632 

8.06xlO“ 10 

55.55 

28.9 

17.8 - 47.8 

3FGL J1225.9+2953 

1.42xl0 -9 

83.76 

47.8 

30.0 - 64.9 

3FGL J1315.7-0732 

8.35xlO“ 10 

54.83 

34.6 

17.8 - 77.7 

3FGL J1543.5-0244 

7.35xlO“ 10 

38.90 

10.0 

< 40.7 

3FGL J1544.6-1125 

l.OlxlO -9 

32.98 

12.0 

< 16.8 

3FGL J1601.9+2306 

8.56xlO“ 10 

46.94 

27.4 

< 187.7 

3FGL J1625.1-0021 

3.6xl0“ 9 

31.84 

35.9 

31.6 - 41.4 

3FGL J1659.0-0142 

9.26xlO“ 10 

23.91 

37.9 

< 273.8 

3FGL J1704.4-0528 

8.89xlO“ 10 

20.80 

35.2 

< 964.7 

3FGL J1720.7+0711 

8.89xlO“ 10 

23.41 

22.9 

< 72.3 

3FGL J1744.1-7619 

3.85xl0“ 9 

-22.47 

28.4 

25.5 - 32.2 

3FGL J1803.3-6706 

7.12xlO“ 10 

-20.37 

45.3 

25.0 - 72.4 

3FGL J1946.4-5403 

1.72xl0 -9 

-29.56 

22.1 

17.2 - 26.9 

3FGL J2112.5-3044 

3.26xl0“ 9 

-42.45 

51.4 

43.7 - 58.3 

3FGL J2134.5-2131 

7.06xlO“ 10 

-45.08 

50.5 

17.5 - 165.5 

3FGL J2212.5+0703 

1.24xl0 -9 

-38.56 

34.0 

21.8 - 51.5 

3FGL J2103.7-1113 

1.091xl0 -9 

-34.42 

31.6 

21.3 - 47.0 

3FGL J2133.0-6433 

8.36xlO“ 10 

-41.27 

26.4 

13.1 - 61.5 


TABLE I. A list of Fermi’s most promising dark matter subhalo candidates. In particular, this table includes all unassociated, 
non-variable 3FGL sources with \b\ > 20°, <f> 7 > 7 x 10“ 10 cm -2 s -1 (> 1 GeV), and that are well-fit (% 2 < 25) by annihilating 
dark matter for at least one value of the dark matter’s mass (assuming annihilations to bb). 



FIG. 5. The distribution in galactic latitude and best-fit dark 
matter mass for those bright (F 1 > 1.0 x 10 -9 cm _2 s _1 , > 

1 GeV), non-variable (variability index < 80), unassociated 
3FGL sources located outside of the Galactic Plane (|6| > 5°). 
The error bars around each point represent the range of dark 
matter masses for which the fit to the spectrum yields y 2 < 
25. In this figure, we have assumed annihilations that proceed 
to bb. 


lations to bb). In our opinion, the sources contained in 
this list merit further investigation. If associated emis¬ 
sion can be detected at other wavelengths, these sources 
could be excluded as subhalo candidates, allowing us to 
derive stronger constraints on the dark matter annihila¬ 
tion cross section (see Sec. VI). Alternatively, the lack of 
counterparts at other wavelengths could strengthen the 
case that one or more of these candidates are, in fact, 
dark matter subhalos. 


V. SPATIAL EXTENSION 

Most astrophysical gamma-ray emitters (pulsars, blaz¬ 
ers, etc.) are effectively point sources, without any po¬ 
tentially observable spatial extension. A fraction of dark 
matter subhalos, on the other hand, could be detectably 
extended, especially those most nearby and large. In 
this section, we study the dark matter subhalo candidate 
sources identified in the previous section in an effort to 
determine whether they exhibit any evidence of spatial 
extension. 

The Fermi Collaboration has previously performed 
a search for spatial extension among their brightest 
sources, using their first two years of data |23^ . In do- 
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FIG. 6. The change in the test statistic (TS) for Fermi’s unassociated, very bright (P 7 > 10 9 cm 2 s 1 , E 1 > 1 GeV), high 
latitude (|6| > 20°) sources when the point-like template is replaced with a spatial template with gaussian extension. 


ing so, they identified 21 spatially extended gamma-ray 
sources, none of which appear likely to be dark matter 
subhalos (17 of these 21 objects lie within 10 degrees of 
the Galactic Plane, and the four others are each associ¬ 
ated with emission from the Large and Small Magellanic 
Clouds, the nearby galaxy Centaurus A, and the Ophi- 
uchus molecular cloud). Given the fact that they did not 
identify other spatially extended sources, we do not ex¬ 
pect any of the subhalo candidates under consideration 
to exhibit extension at a highly statistically significant 


level. Given the larger data set that is currently avail¬ 
able, however, it is possible that spatial extension could 
be detectable in one or more of our subhalo candidates 
at a non-negligible level. 

To test for evidence of spatial extension, we repeat 
the procedure described in the previous section, replac¬ 
ing the point-source template with a gaussian template 
of width a. In Fig. [6j we plot the change in the test 
statistic (TS) (defined as twice the difference in the 
global log-likelihood) of the source when the point-like 
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J2212.5+0703 



FIG. 7. The change in the test statistic (TS) for Fermi’s unas¬ 
sociated source J2212.5+0703 when the point-like template is 
replaced with a spatial template with gaussian extension. 


template is replaced with that of an extended source. 
Of these 18 sources considered (those with |6| > 20° 
and F 7 > 10 -9 cm -2 s _1 , E 1 > 1 GeV), we found 
one for which the TS increased when a spatially ex¬ 
tended template was adopted (see Fig. [6]). This source, 
J2212.5+0703, is one of the promising high-latitude dark 
matter subhalo candidates identified in Sec. [TV| and is 
well fit by m x ~ 21.8 — 51.5 GeV (for annihilations to 
bb); see Table [ij We find that the TS of this source in¬ 
creases by 17.9 when the point-like template is replaced 
by a gaussian template of width 0.15°, corresponding to a 
local statistical significance of 4.2cr (or 3.6cr after account¬ 
ing for a trials factor of 12, corresponding to the number 



FIG. 8. The spectrum of J2212.5+0703, compared to that 
predicted from a 34 GeV dark matter particle annihilating to 
bb (dashed curve). 


of subhalo candidate sources tested for spatial extension). 
This is shown in greater detail in Fig. [7] In Fig. [8j we 
show the gamma-ray spectrum of this source. If the sig¬ 
nificance of this extension were to increase as Fermi con¬ 
tinues to collect more data, it would help to support a 
dark matter subhalo interpretation for this source, over 
that of a pulsar or other point-like object. 

As a control group, we looked for signs of spatial ex¬ 
tension among a subset of Fermi sources that have been 
identified as pulsars (see Fig. |9|. In particular, we con¬ 
sider pulsars with the same range of latitudes (|b| > 20°) 
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FIG. 9. As in Fig. [6] but for Fermi sources that have been identified as pulsars that are located at high latitude (\b\ > 20°) 
with a gamma-ray flux in the range of + 7 = (1.0 — 4.3) x 10 -9 cm -2 s -1 (F? 7 > 1 GeV). No evidence of spatial extension is 
found among this control group of sources. 


and gamma-ray fluxes (F 7 = 1.0 — 4.3 x 10 9 cm 2 
s _1 , E 1 > 1 GeV) as those subhalo candidates shown in 
Fig.© We find no evidence of spatial extension among 
this control group of sources (the most significant prefer¬ 
ence for extension was from J1630.2+3733, for which the 
TS increased by 0.6 when modeled with an extension of 
(7 = 0.1°). 
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VI. CONSTRAINING THE DARK MATTER 
ANNIHILATION CROSS SECTION 

In this section, we use the results of Sec.[lV]to place an 
upper limit on the number of bright, high-latitude dark 
matter subhalos, and use this to derive a constraint on 
the dark matter annihilation cross section. We do this 
through the following procedure. First, for each subhalo 
candidate, we calculate the y 2 of the fit for gamma-ray 
spectrum predicted by a given dark matter model, and 
convert this into a p- value. Twice the sum of the p -values 
for all of the candidate sources under consideration rep¬ 
resents the weighted number of sources, which we plot 
in the left frame of Fig. [lOj As already seen in Fig. [5j 
this result indicates that ^10-15 of these sources could 
plausibly be subhalos of ^20-50 GeV dark matter parti¬ 
cles (they could also be pulsars, however). Dark matter 
particles with heavier masses (> 100 GeV), in contrast, 
do not provide a good-fit to any of the subhalo candidate 
sources. 

In the right frame of Fig. [lOj we plot as a solid line the 
95% upper limit on the dark matter annihilation cross 
section derived from this source population. This result 
is obtained by calculating the Poisson errors around the 
weighted number of sources shown in the left frame. The 
dashed line represents the constraint that would have 
been derived if zero subhalo candidates had been ob¬ 
served. For high values of the dark matter mass, the 
weighted number of sources is only slightly greater than 
zero, and these two lines are almost identical. For lower 
masses, in contrast, many subhalo candidates exist and 


the resulting constraint is somewhat weaker. 

In the left frame of Fig. El we show this constraint as 
derived for a variety of dark matter annihilation chan¬ 
nels (t + t - , FF + VF _ , ZZ, 65, cc). We do not consider 
annihilations to e + e - , as these final states are 

more strongly constrained by cosmic ray positron mea¬ 
surements |24 than by gamma-ray observations. In the 
right frame of this figure, we compare the constraints 
presented in this paper to those derived from gamma-ray 
observations of dwarf spheroidal galaxies m (see also, 
Ref. [13.:), the Galactic Center (assuming an NFW halo 
profile) [14], the extragalactic gamma-ray background 
(assuming a reference substructure model, and taking the 
less conservative “sensitivity reach” approach) [T5] (see 
also, Ref. csi), and the Large Magellanic Cloud [25]. 
The constraints presented here are comparably stringent 
to those derived from these other observations, and are 
complementary due to the differing systematic uncertain¬ 
ties involved. 

We also show in the right frame of Fig. [Tl] the re¬ 
gion of parameter space that is able to account for the 
gamma-ray excess observed from the region surrounding 
the Galactic Center [26b32l , as presented in Refs. G3I34]. 
This region assumes a local dark matter density of 0.4 
GeV/cm 3 and a scale radius for the Milky Way’s halo 
profile of 20 kpc. Varying these quantities within their 
current uncertainties can allow this region to shift up¬ 
ward or downward by an additional factor of a few. It 
is interesting to note that if the Galactic Center excess 
is generated by annihilating dark matter, then we expect 
that a few of our dark matter subhalo candidate sources 








13 




FIG. 10. Left: The weighted number of unassociated, non-variable (variability index < 80), bright (F > 7 X 10“ 10 cm -2 s -1 , 
E 1 > 1 GeV), high-latitude (\b\ > 20°) 3FGL sources as a function of dark matter mass (assuming annihilations to bb ), defined 
as twice the sum of the p -values for all of the candidate sources under consideration. Right: The 95% upper limit on the dark 
matter annihilation cross section derived from this source population. The dashed line represents the constraint that would 
have been derived if zero subhalo candidates had been observed. For high values of the dark matter mass, the weighted number 
of sources is only slightly greater than zero, and these two lines are almost identical. For lower masses, in contrast, many 
subhalo candidates exist and the resulting constraint is somewhat weaker. 




FIG. 11. Left: The 95% upper limit on the dark matter annihilation cross section, for a variety of annihilation channels. Right: 
A comparison of the limit presented here to those derived from gamma-ray observations of dwarf spheroidal galaxies |12| . the 
Galactic Center m, the extragalactic gamma-ray background m, and the Large Magellanic Cloud EH- Also shown is the 
region of parameter space that is able to account for the gamma-ray excess observed from the region surrounding the Galactic 
Center j26ll32] . as presented in Refs. [33][34] (assuming a local dark matter density of 0.4 GeV/cm 3 and a scale radius for the 
Milky Way’s halo profile of 20 kpc). 
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should, in fact, be dark matter subhalos. 


VII. SUMMARY AND PROSPECTS 

The dark matter halo of the Milky Way is predicted to 
contain a very large number of smaller subhalos, and if 
the dark matter consists of weak-scale annihilating parti¬ 
cles, the most massive and nearby of these subhalos could 
appear as gamma-ray sources without counterparts at 
other wavelengths. In this paper, we have studied the 992 
unassociated gamma-ray sources contained in the Fermi 
Collaboration’s Third Source Catalog (3FGL), in an at¬ 
tempt to constrain the number of these sources that could 
be dark matter subhalos. From this information, we have 
derived constraints on the dark matter annihilation cross 
section that are comparably stringent to those based on 
gamma-ray observations of dwarf spheroidal galaxies, the 
Galactic Center, and the extragalactic gamma-ray back¬ 
ground. 

We also present in this paper a list of 24 sources 
which we consider to be promising dark matter sub¬ 
halo candidates. These sources are each quite bright 
(F 7 > 7 x 10 -10 cm -2 s -1 , Ery > l GeV), are located 
far away from the Galactic Plane (\b\ > 20°), show no 
signs of variability (variability index < 80), and exhibit 
a spectral shape that is compatible with that predicted 
from annihilating dark matter particles. The majority of 
these subhalo candidate sources are best-fit by dark mat¬ 
ter particles with masses in the range of ~20-50 GeV (as¬ 
suming annihilations to bb ), similar to that favored to ex¬ 
plain the previously reported Galactic Center gamma-ray 
excess [26d34] . Gamma-ray pulsars are known to exhibit 
similar spectra, however, making it difficult to determine 
whether these sources are dark matter subhalos or radio¬ 


faint gamma-ray pulsars. We also report 3.6cr (global) ev¬ 
idence of spatial extension from the bright, high-latitude 
subhalo candidate source 3FGL J2212.5+0703. If the sig¬ 
nificance of this extension continues to increase as Fermi 
collects more data, and no multi-wavelength counter¬ 
parts are identified, it would strengthen the case that 
this source could be a dark matter subhalo, rather than 
a pulsar or other point-like gamma-ray source. 

Of the 117 high-confidence gamma-ray pulsars de¬ 
scribed in the Second Fermi-LAT Pulsar Catalog [35] . 
36 have pulsations that were first identified in gamma- 
rays (as opposed to in radio or other wavelengths). Each 
of these sources has been studied with deep radio ob¬ 
servations f36H38] . yielding several detections 
The fact that dozens of gamma-ray pulsars have not yet 
been detected at other wavelengths, however, leads us to 
expect that some of Fermi’s unassociated sources could 
be radio-quiet pulsars, whose gamma-ray pulsations have 
thus far gone undetected (perhaps due to weak pulsation 
and/or broad pulse profiles). Deeper multi-wavelength 
observations and further searches for gamma-ray pulsa¬ 
tions will be essential if we are to identify the nature of 
these intriguing sources. 
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